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1. Algebraic Preliminaries 

The goal of the present paper is primarily to exhibit the effectiveness of using 
Lie algebras to compute explicit perturbation eigenvalues for quantum anharmonic 
oscillators in one dimension. There arc however, several goals secondary in stature, 
but which merit discussion. The first of these is to enable the reader to work with 
Weyl algebras in the abstract. Two presentations of it arise readily in quantum 
physics. In particular the first is the algebra of position and momentum operators 
in nonrelativistic mechanics wherein [x,p\ — ih. The second is the algebra of ladder 
operators [a, = 1. It is the second presentation with which we will be primarily 
concerned in this paper. Of course, the first presentation may be made to look like 
the second by considering not p, but instead the simple derivative ^ whereby one 
has = 1. 

1.1. Normal Ordering and Weyl Binomial Coefficients. For any abstract 
Weyl algebra determined by two elements A and B obeying [A, B] = 1, an ordering 
of a polynomial in A and B will be said to be normally ordered if all powers of B 
appear to the left of powers of A. For example A^B^ is not normally ordered, but 
B'^A^ is. In our case a) will always be placed to the left of a. As it is well known in 
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elementary quantum mechanics one may move back and forth between presentations 
of problems in position-momentum coordinates and annihilator-creator coordinates 
with the following equivalences 

n 1 1^ a + a'f 

(1.1.1) x = 



V2 

a 

P = 



a — 



iV2 

Since this paper is concerned with anharmonic oscillators we will be concerned with 
in the potential. Thus, we need an efficient way of normally ordering (a + a^)". 

Lemma 1. Let A and B determine a Weyl algebra so that [A,B] = 1. The normal 
ordering of {A + B)"' is given by 

n min{n,n— m} , >^ 

(1.1.2) {A+By = Y^ \ TO r 

rTi=0 ;c=0 ^ 



where 
(1.1.3) 



2^k\{m - k)\{n - m - k)\ 

is the Weyl binomial coefficient. 

The proof of this lemma involves nothing more than counting commutations. 

Example 2. We will use the fourth order relation explicitly later, so here is an 
example of how the Weyl coefficients factor in. 

(a + a'^f = a^f^ + Aa'f^a + Ga^'^a'^ + 4a^a^ + + Qa''^ + Ua^^a + 6a^ + 3. 

Remark 3. Notice that 



r n 1 




r n \ 


I J 




^ n — m j 



If one wishes to attempt calculations within a Weyl algebra it may be useful to 

compute with abstract elements A, B first and then plug into a specific situation 
one has in mind. One other useful tip is that if one has Weyl variables A, B then 
it can be convenient to consider representing the algebra as B or ^, — This 
becomes consistent with the first presentation considered. For example, one peculiar 
formula which is arguably easier to compute with abstract Weyl variables is 

(1.1.4) f,-^^f{x)=f{fix). 

1.2. Baker-Campbell-HausdorfF and the Hadamard Lemma. We will be 
concerned throughout much of this paper with exponentiating noncommuting vari- 
ables. We run into a stopping block in trying to compute the exponentials explicitly. 
The main issue is that for noncommuting variables X, Y we see 

The Baker-Campbell-Hausdorjf formula is the solution to Z = log(e"^e^). The 
explicit solution is formally given as symmetric sums and differences of nested 
commutators in X and Y. One may find this expression in nearly any textbook on 
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advanced quantum mechanics. We shall not be concerned, however, with isolated 
exponentials, but rather expressions of the form 

Using elementary combinatorics and the Baker-Campbell-Hausdorff formula one 
can arrive at the Hadamard lemma. 

Lemma 4. Let X, Y be noncommuting variables then one has 

(1.2.1) e^Fe-^ =Y +[X,Y] + [X,Y]] + [X, [X,Y]]] + ■■■ 

If we allow the notation [X^'^\Y\ = [X, [. . . , [X, Y]\\ then we may write more suc- 
cinctly 

(1.2.2) e^Fe-^ =^^[x('=),r]. 

fe=o 

1.3. The Formula Often Desired and Rarely Known. One final assertion 
about Weyl variables in the algebraic preliminaries must be the formula 

min{n,m} , \ / \ 

(1.3.1) [A-,B^] = B'^-^A^-K 

This formula is often loft as an exercise in quantum mechanics texts and some- 
times in homological algebra, but rarely is it completed. One might jokingly say 
it is similar to the snake lemma in that no one knows if it's really true since the 
only persons who have ever proven it are graduate students. All kidding aside, 
this is indeed the correct formula for normally ordering variables obeying the Weyl 
relation. 

2. The Method of [JA] 

The main impetus for this research comes from the paper [J A] . The goal of this 
section is to explicate in a reasonably clear manner the content of that paper. 

The premise upon which [JA] begins is the idea that we can create a new Lie 
algebra from simply taking commutators of the unpertTirbcd Hamiltonian and 
the new anharmonic Hamiltonian ff„ . For the sake of mathematical simplicity the 
Hamiltonians in question are given as essentially unitless operators: 

(2.0.2) Ho= i(p2+a;2)= ata+^ 

As one may infer, we have made the following assumptions and simplifications: 
(1) h = LO = m = 1, 

V2 ~ V2 ' 

Let us now give the formulation of the Lie algebras. 
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Definition 5. The Lie Algebra An — {Lm}rnei is generated by the elements 

Li = Hq, L2 = Hn 

and other satisfying 

(2.0.3) [Li, Lj] = CijmLjn- 

for some structure constants Cijm € C. Furthermore, this Lie algebra should be 
closed under commutators up to order A*^. In other words no Lm should be of the 
form A'^+^(a^*a'' — a^*a*) for any s,t. That is, formally we require A*^"*"^ = within 
the Lie algebra. 

In the case of this paper we will consider An ^ unless otherwise explicitly stated. 
In fact, [JA] only considers Lie algebras up to order one in A with the exceptions 
of n = 1 and n = 2 because these determine harmonic oscillators and their solu- 
tions are already known. We deal with the special technique for solving harmonic 
oscillators in the appendix. 

Once the algebra An^ is determined we proceed in the following way. Suppose 

i 

(2.0.4) [Li,X2] = ^ci2feLfe 

k=3 

where each L). is of the form 

(2.0.5) A(at™a^ - at^a™). 

The symmetry of these Lk is important and comes back in an important way 
due to the normal ordering procedures we have adopted. We will see this explicitly 
in the computations. 

We then construct a unitary element of the associated Lie group by 

3 

(2.0.6) U = exp{J2^kLk). 

fe=3 

This says that the only Lk allowed in our unitary are those arising directly from 
the commutator [Li, L2]. The ak are real constants which we will tune as necessary. 
Once we produce such a unitary we make a transformation from Hq to Hn by 

(2.0.7) U^HoU = F„ - A„. 

In each case A„ is an operator which simply controls the pertubations of eigenvalues. 
Furthermore, by the clever choice of U we will have [C/, A„] = + O(A^). Due to 
the Hadamard lemma we can produce A„ by computing simple commutators. 

At this stage one may write the new eigenvectors as W\j) where |j) are the 
eigenvectors for the harmonic Hamiltonian with eigenvalues i + 5- Therefore, up 
to order A^ our equation now reads 

(2.0.8) HnU^\j) = {U^HoU + An)U^\j) 

= t/tffo|j)+A„C/t|j) 

= U\j+^-)\j} + U^Ar,\j} 

= {j + l + K)u^\j)- 
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In essense, depending on the form of A„, we will be able to road off the first 
order perturbation eigenvalues (A„) of Hn with relative ease. 

A natural question arises as to when we can solve this system explicitly. [JA] 
makes a passing statement which we will now state as a formal theorem. 

Theorem 6. If the Lie Algebra An is closed (in all orders of X) then we can solve 
the n*'' order anharmonic oscillator in closed form. 

Proof. Let An = {Lk}k=i be a closed Lie algebra corresponding to the Hamiltonian 
Hn- Then consider the general Lie group element given by 

N 

U = exp(^ akLk) =: exp(L). 

k=l 

From the Hadamard lemma we obtain 
(2.0.9) U^HoU = Y,^[L^''\Ho]. 

k=l 

Since An is closed, the commutators [L^^^ , Hq] either vanish or give Lie algebra 

elements with some periodicity. In this way we can formally sum them in a power 
series. Setting our parameters to appropriate values we obtain 

U^HqU — H„ + perturbations. 

□ 

3. Lie Algebras up to Order One in A 

We begin the computation of the Lie algebras by giving an important commu- 
tator relation. 

(3.0.10) [a^a, a^'^a^ ± a^^a''] = (k - e)ia^''a^ T a^^a'') 

This equation paired with the symmetry of Weyl binomial coefficients points us 
to assigning a^^a"^ — a^'^a"' as our Lie algebra elements. Let us begin by computing 
[Hq, Hn]. 

[Ho,Hn] = [Ho,Ho + ^{a^ +ar] = 

v2 

[Ho,-^{a^+ar] = -|.[ata,(at+a)"] 

A 



* k.m 



a 

k ' 



(3.0.11) = }^(2"^-«)(« 



fm— /c^n— m— ^fn— m— fe^r) 



We will throw away the multiplicative constants in favor of rescaling them by ak 
in our general Lie group element. Therefore, the first batch of elements revealed 
to us are those of the form A(a^^a'" — a^^a^). Once we realize these, we begin 
commuting again with Hq to find more elements of the form A(a^^a™ + a^™a^). In 
order to close An'' we also need to add a central element / to our Lie algebra. All 
other commutators will involve terms with A'^ and therefore we disregard them in 
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Remark 7. Notice that no elements of the form a^"a^ appear anywhere. This 
is because they can be written in terms of the number operator N = a)a which 
commutes with _ffo- 

Given a general Hamiltonian with special exceptions, by simple combinato- 
rial formulae one infers the number of generators for A^n'' by 

(3.0.12) |4fc| = (fc + l)fc + 3 

l4lVil = (fc + l)A; + 3 

Example 8. Let's take a quick look at ^g^' . Of course, we let Li = Hq and L2 = Hq. 
By our computation we know the generators arising from [Li, L2] are as follows: 

A(at6-a6), X{a^^a-a^a^), A(atV - at V), 
A(at4-a4), A(at3a-ata3), A(at2 - a^). 

Furthermore, commuting these with Li we obtain simlarly symmetric elements 
with plus signs. Finally we add in I to account for commuting elements. Notice 

if we commute any other elements we obtain an element in which we have 

formally disallowed for now. Therefore |^g^^ | = 15 = (3 + 1)3 + 3 as previously 
stated. 

4. Explicit Computations 

In this section we will derive the first order perturbation for all anharmonic 
oscillators with Hamiltonians of the form 

There are two distinct cases for computing first order perturbations; odd and even. 
We will treat the odd case first. 

For the sake of uniformity in our calculations we will consider Hamiltonians of 
the form H2k-i and H2k- 

4.1. Odd Powered Potentials. From our earlier computations of [Hq, H2k-i] 
and our prescribed form of U wc have 

m-l A- 

(4.1.1) U = exp(A J2 E «-,^a^"""'"'«' - at^a^""'"'). 

e=o m=l 

By requiring am,e € M we obtain W = and we may now apply the Hadamard 
lemma. Letting 

rn— 1 k 

e=o m=l 

we have: 



(4.1.2) 



U^HoU = Ho + [-X, Ho] + h-X, [-X, Hq]] + ■ 
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We notice immediately that X contains a multiplicative factor of A and since we 
have = we may ignore all terms past [—X, Ho] . 

Using previous computations and elementary properties of derivations we have 

m— 1 k 

(4.1.3) [Hq, X\ = \^Y^ a^,e{2m - 1 - 2£)(a^2m-i-£^£ ^ a^^a^"-!-^). 

e=Q m=l 

If we recognize that 

(nJ-nU'^k-l 1 2fe-lmin{2fe-l,2fe-l-m} 



X = 



m=0 j=0 ^0 



we see that in order to produce H2k-i = -ff o + Aa;^*' ^ we need to set 

(4.1.4) a.,(2m-2£-l) = -^{ }^_^. 

Notice what wc have done. We have transformed Hq into H2k-i+0{X^). Hence, 
there is no perturbation term up to first order. 

Example 9. Let us compute the example of Hi = Hq + Xx explicitly. We already 
know that this is a shifted harmonic oscillator, where a simple change of variables 

reveals the energy eigenvalues arc n + ^ — 

In our case the Lie algebra simplifies slightly and is fully closed as 

Ho,Hi,X{a^ - a), I. 
Our appropriate unitary transformation U is therefore given as 

U = exp(Q;A(a^ ~ o-))- 
We can compute this even more explicitly than before given that [a^ — a,a^ = 
[a' — a,a] = —1. 
In this case we know 

(4.1.5) [A,B]=I3gC =^ [A, e^] = pe^ 

for any abstract operators A, B. 

Therefore [a, U] = aXU and [a^,U] = aXU 

U^HoU = U^a^aU+^ 

= U^a^U{aX + a) + ^ 

= U^U{a^ +aX){a + aX) + ^ 

= a^a+^+ aX{a^ + a) + a^X^. 

Setting a = we derive 

(4.1.6) U'^HoU = Hi + —. 
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This is exactly the result we previously knew. Notice however, there is no A term. 
The first perturbation term is order A^. 

Example 10. To better see the odd powered result more explicitly, let us compute 

the result for H^. 

Our Lie algebra up to order one is given by 

A(at5-a5), X{a^^a-aU^), X{a^^a^ - a^^a^), 
A(at3-a3), A(at2a-ata2), A(at - a), 

A(ot3 + a3), A(at2a + ata2), X{a^ + a), 
Ho, H5, I. 

Therefore U is 

U = exp(A(a3,o(a^^ - a^) + 0:3,1(0^^0 - a^a^) H h ai,o(a^ - a))). 

Our Hamiltonian transforms as 
U^HoU = Ho + X[Ho, a3,o(a^' - a^) + a3,i{a^*a - a^a^) + ■■■+ ai,o(at - a)] + OiX") 
= Ho + A(5Q!3,o(a^^ + a^) + 3a3,i{a^'^a + aU'^) + ■■■ + ai,o(a^ + a)). 

Setting 



"3,0 


= 2-5/2/5, 


Q!3,l 


= 2-5/25/3, 


as, 2 


= 2-5/210, 


^2,0 


= 2-5/210/3, 


a2,i 


= 2-5/230, 


ai,o 


= 2-5/^15, 



we obtain 



(4.1.7) U^HoU = Ho + Xx^ + 0{X^) = Hr, + O(A^). 

4.2. Even Powered Potentials. It is the goal of this section to show the first 
order perturbation energies for oscillators corresponding to i?2fc are 



(4.2.1) 




This result is obtained rather easily utilizing the technology we have developed 
for odd powered potentials. We only need to realize that the operators which are 
not Lie algebra elements are of the form a^"a". This comes from the required 
symmetry of our generators. Hence our unitary U will appear exactly as in the odd 
powered case and the Hadamard lemma yields: 

(4.2.2) U^HoU = Ho + Xx^'^ - ^ ^ | ^ M ^ti^j + o(A2). 

The only thing left to pretty up our example is changing a^'^a'' into an expression 
of number operators. 
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Recall N = a^a has nonncgative integer eigenvalues given hy N\n) = n\n). In 
this way any expression f{N) in our perturbation expansion will give eigenvalues 
f{n) by the functional calculus. 

Proposition 11. 

(4.2.3) atfea'^ = k\(^^ 

Proof. We refer back to our commutation relation from §1.3. Thus we 

have 

= at(aat'=-i-(fc-l)at'=-2)a'^-i 

= (ata)atfe-ia'=-i_(fc_i)at'=-ia'=-i 

= {N -{k-l))a^''-'a''-\ 



N 
k 



Repeating this we see 

(4.2.4) a^'^a'^ = N{N - I) ■ ■ ■ {N - {k - I)) = k\ 

It is merely a matter of rearranging terms to see 

(4.2.5) U^H,U = H,,-^^j2j\{ ] { ^ )+0(A^). 
Example 12. Let us look briefly at the Hamiltonian 



□ 



F4 = a'fa + i + ^(a^ + a) 



4 



This is the famous quartic which has received much attention in texts and papers. 
We can check our results against those of standard perturbation theory. 
Our Lie algebra ^4^^ is given by 

A(at4_a4)^ A(at3a_ata3)^ A(at2 - a^), 
A(at4 + a4), A(at3a + ata3), X{a^^+a^), 
Ho, Hi, I. 

Our unitary is given by 

U = exp(A(a2,o(a^^ — a^) + Q!2,i(a^^a — a^a^) + aifi{a' — a))). 
Prom here we must simply crank the handle for our machine and we realize 
U^HoU = i/o + [^^o,A(a2,o(a^^-a4) + a2,i(a^^a"a^a3)+ai,o(a^-a))] 
(4.2.6) = ifo + A(4a2,o(a^^ + o^) + 2Q2,i(a^^o + ota3) + ai,o(at +0)). 
Setting 

"2,0 = 1/16 
"2,1 = 1/2 
ai,o = 3/4 

we arrive at 
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(4.2.7) U^HqU = i?4- ^(6a^^a2 + 12a^a + 3) 

= ff,-^(iV(Ar + l))-^. 
If we look to the ground state we see that 

which agrees with the standard perturbation theory. 

In particular the perturbed ground state of the anharmonic oscillator correspond- 
ing to is given by 

5. Extending the Method 

The second subsidiary goal of this paper is to show several extensions to this 
method and invite research into even more applications of Lie algebras into physics. 

5.1. Simple One Dimensional Corollaries. Now that we have given explicit 
formulae for computing perturbation eigenvalues for potentials of the form Aa;" 
we can extend by linearity (up to order one) and immediately recover eigenvalues 
for polynomial potentials. In fact we can extend this further to convergent power 
series. 



Example 13. Let us consider 



Of course this can be rewritten 

k 

If we notice that only the even powered potentials contribute perturbations up 
to first order then we will also compute perturbations for H = Hq+X cosh(x) as well. 

Let us compute only the ground state energy. We have 

= 5 + A(E = i + A(i: = 1 + ,<„p,v4). 

k ^ ' k 

Moreoever, we can add any number of perturbation parameters and solve the 
system accordingly. In particular we can essentially read off first order perturbations 
for Hamiltonians of the form 

n 
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5.2. Simple N Dimensional Corollaries. In extending this method, it is natural 
to ask whether one can tackle higher dimensional systems with a similar approach. 
In our case, we certainly can attack higher dimensional problems similarly, but the 
construction of the Lie algebra is different. For the simple A'' dimensional corollar- 
ies, we will assume our oscillator potential is not coupled (i.e. No terms of the form 
Xx^y'') appear. For the sake of simplicity let us go through the construction of the 
Lie algebras for a two dimensional oscillator. 

Consider ^ ^ 

Hn,m = alax + 2 + '^i-'^" + "l^-v + 2 ^ -^22/™- 
We will take four elements as given in our Lie algebra 

t 1 t 1 
Hofi = flTfl^ + - + alay + -, 

Hn,m = iIo,o + Aia;" + A22/". 

In this way we will set up our Lie algebra as two independent oscillator Lie 
algebras and solve our problems from before. Consider for example 

-ffi,4 = Hofi + Aix + A2t/''. 

Our Lie algebra A'^iJ^^ will have the following elements 

Ho,o, Hifi, Ai(o|. — Ox), 

HoAi Hi4, I, 

A2(4'' + a*), X2{al^ay + alal), X2{al^ + al). 

Our unitary takes the form 

U = cxp(aAi(4 - a^) + f3iX2{aP - a^) + /^aAal^^Oj, - 4^) + AiAal^^ - 4)). 

Now we use the Hadamard lemma again, but taking advantage of the relations 

(5.2.1) [a],ak] = Sjk 

we can completely separate x variables from y variables and our calculation plays 
out exactly as before. 

For the Hamiltonian Hi ^4 our perturbed ground state is 

Eo = l + 0{Xl) + l + ^-^ + 0{Xl). 
Now we can use all the simple one dimensional corollaries in turn as well. 

5.3. Higher Order Perturbations. Since perturbation theory is meant to com- 
pute more than first order terms we seek to use this Lie algebraic method to compute 
higher order terms. Certainly, one can see that using the transformations U we have 
set up thus far will produce higher order terms. One can see this if we set 

U = exp(AL). 

Our transformation becomes 
(5.3.1) U'fHoU = Ho + X[Ho,L] - ^[L, [Ho,L]] . . . 
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This approach, however, changes our Hamiltonian funadmentally. In fact, we end 
up not solving any problems, but instead creating more. A quick trial calculation 
with any Hamiltonian carrying term or higher will reveal that we cannot cancel 
certain terms arising from [L, [I/,i?o]]- To remove this difHculty we must expand 
our Lie algebra to include terms carrying A'^ for whichever k we should choose. It 
is therefore convenient to write our unitary transformation as 

fc 

(5.3.2) [/ = exp(^A^L(,)), 

where iy) are Lie algebra elements arising from order comutators. 

Example 14. Let us return briefly to the quartic oscillator and calculate its second 
order perturbation. Since it is well studied we may verify our results easily. 

Computing commutators and commutators of commutators one will arrive at 
the following Lie algebra up to order 2. 

Ho, Hi, I, 

A2(at6±a6), A2(at4a2±atV), 

Knowing the form of our necessary first order transformation, we add four terms 
to the exponential by 

(5.3.3) U = exp(A(47(a^*-a^) + ^(a^^a-a^a^) + ?(a^'-a')) 

16 2 4 

+A2/3i(at6 - a^) + X^p^ia^^a" - a^^a^) 
+X^P^{a)^a - a+a^) + X^ pi{a)'' - a")). 
U = exp(AL(i) + A2L(2)). 

For simplicity, let us compute only the ground state energy given by U^HqU. 
We have 

U'fHoU = exp(-AL(i) - X^L^2))Ho exp(AL(i) + X^L^2)) 

= Ho + X[Ho, + A^ [Ho, L(2)] - y [i(i) , [Ho, L(i)]] + 0(A3) 

Prom here it is a matter of computing commutators and adjusting /3i, /?2j /^a, (^4 
to cancel higher order terms not given as functions of number operators. 

When we compute the ground state energy we are concerned only with constant 
terms. Therefore, looking to our commutators we have the A^ term 

[^(at4-a4) + l(at3a-ata3)+3^«t2_^2^^1(^t4+^4^ + (^t3^+^t^3)+3^^t2^^2)] 

Expanding this we are left with two terms giving constants 

j-[a^^-a\a^^ + a^] and ^[a^^ - a\ a^^ + a% 
64 8 

Our constant terms turn out to be — and "^g^'-*^ yielding 
Finally, our ground state energy up to second order will be given by 
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(5.3.4) {U^HoU)U^\0) = {H,--- _-)[/t|0), 
yielding 

(5.3.5) Eo = l + ^-^-^+0{X^). 

Indeed, this ground state energy agrees with the standard perturbation theory. 



For the interested reader, the correct /3 parameter values are 

1 

48' " 



/3i - ^> - -f|, 
-9 



and the second order equation appears as 

(5^«) 5U^(«)+H^(-^)+36A=« + a^+0(A3, 

6. Discussion 

One important problem among many arising from this paper and which we have 
yet neglected to mention is the representation theory the Lie algebras. For example 
if we are dealing with the cases Xx,Xx^ in the one dimensional case or any qua- 
dratic term in higher dimensional cases, wc have a closed Lie algebra. This may 
not be terribly surprising as a closed Lie algebra offers an exact solution, and these 
particular potentials are simply shifted or coupled harmonic oscillators. However, 
the Lie algebras wc have constructed all contain central elements. In the cases of 
Ai and A2 we have 4 dimensional closed Lie algebras with center. Representation 
theory tells us that these arc isomorphic to 9(2. It remains to be seen exactly the 
relationship between these Lie algebras and the symmetries they describe. In this 
paper we have only used them to calculate perturbed energy levels. It is entirely 
possible there is a simpler approach to this problem from an entirely representation 
theoretic standpoint. 

As it stands, we have; given explicit constructions for Lie algebras up to any 
order and the method by which we may construct a unitary operator to make the 
transformation 

Hq Hn + perturbations up to (9(A'^). 

By taking advantage of our symmetric construction of these Lie algebras, the 
Hadamard lemma, and several formulae concerning abstract Weyl algebras we have 
managed to give eigenvalues in agreement with standard methods. 

Another issue which we have neglected to resolve it how to deal with coupled 
oscillators in general. In the appendix we briefly mention the way to deal with 
potentials of the form Xxy. The Lie algebra computations for coupling terms of the 
form Xx"y"^ are more taxing and trickier. This method, however, should be able 
to deal with situation, but some coordinate change may be required first. 
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It is the hope of the author that anharmonic oscillators are simply a useful class 
of examples for the propitiation of this method. Furthermore, it is hoped that 
this method will help to give rise to additional representation theoretic methods in 
physics. 

Appendix: Dealing with Harmonic Oscillators 

Two important cases we haven't touched upon are those of actual harmonic oscil- 
lators in one and multiple dimensions. Consider for example the two Hamiltonians 

(6.0.7) = a^a+\ + X{^f 

He = al.ax+a\iay + l + ^{a^^ + ax){a\ + ay). 

These correspond to the shifted frequency oscillator in one dimension with new 

frequency \/l + 2\ and a coupled oscillator in two dimensions with quadratic cou- 
pling term. These cases have been well studied and so we have neglected them thus 
far. However, the techniques to compute the perturbations are special because 
these Hamiltonians along with and iJo.o produce closed Lie algebras. By our 
theorem earlier we know that we can solve these exactly and not concern ourselves 
with k*^ order perturbations. 

The main technique we employ is to transform our ladder operators via the so 
called Bogoliubov tranforms. In one dimension we have 

(6.0.8) 6^ = U'^a^U = aa^ +Ta 

b = aU = aa + TO^ . 

In this way we produce the new Hamiltonian 

y/lT2X{b^b + 1) = ata + i + ^(at + af. 

This algebra to move from a, r to this clean form of the new Hamiltonian is tedious 
to be sure. The interested reader should confer with [JA] or email the author for a 
small set of notes. 

A similar technique can be used for the quadratic coupling, but the transforma- 
tion must take into account much more coupling. Our transformation should look 
something like 



(6.0.9) 
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This is simply a coordinate change which decouples the coordinate variables. 
The matrix, however, will take a very special form so that [bi, bj] = Sij as did the 
initial coordinates. 

Remark 15. Notice here that we can couple our ladder operators in many more 
ways. For example we can tackle problems such as dynamic coupling 

H = x''-dl + y^-d^ + \d,dy, 
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or oscillators in a magnetic field 



H ^x' -dt + y' -dt,+ \{yd^ - xd. 



Vh 

So long as our coupling term contains terms of order two or less in each of the 
ladder operators we can tackle thess problems with a simple coordinate change. 
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